Arguments towards a c-theorem from branch-point twist fields 

Olalla A. Castro- Alvaredo*, Benjamin Doyon° and Emanuele Levi* 

** Centre for Mathematical Science, City University London, 

Northampton Square, London EC IV OHB, UK 

° Department of Mathematics, King's College London, 

Strand, London WC2R 2LS, UK 



A fundamental quantity in 1+1 dimensional quantum field theories is Zamolodchikov's c-function. 
A function of a renormalization group distance parameter r that interpolates between UV and 
. IR fixed points, its value is usually interpreted as a measure of the number of degrees of freedom 

^vq I of a model at a particular energy scale. The c-theorem establishes that c{r) is a monotoni- 

cally decreasing function of r and that its derivative c{r) oc — r^(0(r)0(O)) may only vanish at 



O 



quantum critical points (r = and r = oo). At those points c{r) becomes the central charge 
of the conformal field theory which describes the critical point. In this letter we argue that a 
^O ■ different function proposed by Calabrese and Cardy, defined in terms of the two-point function 

(0(r)7'(O)) which involves the branch point twist field T and the trace of the stress-energy 
tensor 0, has exactly the same qualitative features as c{r). 
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1 Introduction 

Branch point twist fields are on the one hand twist fields, that is they are associated to an 
internal symmetry of the theory under consideration, and on the other hand they are related 
to branch points of multi-sheeted Riemann surfaces. The idea of quantum fields associated to 
branch points of Riemann surfaces appeared in [1], where their scaling dimension was evaluated 
in conformal field theory (CFT) (see also [2] for an earlier work concerned with similar ideas). 
Their general description as twist fields associated to a symmetry was given in [3], where they 
were studied in integrable massive quantum field theory (QFT). This description is independent 
of integrability, and it was first used in massive QFT outside of integrability in [3] . 

Branch point twist fields arise quite naturally in the context of the computation of the bi- 
partite entanglement entropy [H [3] . In that context it is convenient to study "replica" versions 
of QFTs consisting of a number n of non-interacting copies of the original model. A twist field 
T is then associated to the Z„ symmetry of the extended model under exchange of the copies. 
More precisely, the field's action is characterized by 

^^{y)T{x) = Tix)^^+liy) x^ > y\ 

^iiy)T{x) = T{x)^^{y) x^<y\ (1.1) 

where ^i{y) represents any field of the theory in copy i, and i = 1, . . . , n with the identification 
n + i = i. With the characterization above the twist field is a local field of the replica theory, 
because the energy density of the replica theory is invariant under permutation of the copies. 
In addition, the usual interpretation of QFTs as perturbed conformal field theories implies that 
any local field of the QFT, including the twist field, has a counterpart in the ultraviolet CFT. 
Branch point twist fields have the further requirement of being spinless and that this counterpart 
be a conformal primary with the smallest scaling dimension. Then, CFT arguments show that 
branch point twist fields have conformal dimensions (A7-, A7-) 

A^ = Ar=^(n-i), (1.2) 

which are functions of the central charge c and the number of copies n [H [2] . 

Since the twist field is associated to a primary field of the underlying CFT, the A-sum rule 
[5] can be employed. This is a rule which provides a means to recover the value ()1.2|) from 
correlation functions in the massive theory. More precisely, it is the statement that the function 
A(r) whose derivative is 

dHr) ^ r((9(r)r(0))(R2)„-(e)(R2)„(r)(R2)J ^^ ^^ 

dr 2(r)(R2),^ 

with A(oo) = 0, satisfies A(0) = A7-. The subindices (M^)„ mean that all correlation functions 
must be considered on n copies of M^, where the twist field is defined. 

The geometrical meaning of the ratio above becomes clear when twist fields are interpreted 
in terms of branch points: 

(e(r)7-(0))(R.)„ _ ^^ ^ ^ ((G(r))^. - (e)R.) , (1.4) 

W/(M2)„ 

where {Q{r))M" is the expectation value of the trace of the stress-energy tensor on an n-sheeted 
Riemann surface with a branch point at the origin and a branch cut in M"*". The various Riemann 
sheets are sequentially joined to each other along the branch cut, as described in detail in [3]. 



In (|1.4p . we used the following steps: 1) in the n-copy model, Q is the sum of the stress-energy 
tensor traces on the various copies, (©(r) • • • )(k2)^ = Yll=i{®j (.''') ' ' ' )(k2)„ and the copies are 
independent, hence (0)(ir2)^ = n{Q)^2; 2) with the twist field insertion, the point (r, 0) on copy 

j is mapped to the point (^^-,0) on sheet j of the Riemann surface — ^~^Y) ^^^ ~ (0(^j))xS 

3) the trace of the stress-energy tensor is spinless, and the theory on the Riemann surface is 
rotation invariant (w.r.t. the origin), so that {Q{rj))M" is independent of j, giving (|1.4p . 

Integrating ()1.3p between and oo we find A(oo) — A(0) = — A7-, a statement that has been 
exploited for example in [1] as an alternative means to deduce the short distance behaviour of 
the entropy in CFT. Integrating (jl.Sp instead between some value tq and 00 we obtain a function 
— A(ro) which "measures" the scaling dimension of the twist field at all energy scales between 
the UV and IR fixed point along the renormalisation group (RG) flow. 

The derivation of (|1.3p carried out in [5] relies heavily on the conservation law of the stress- 
energy tensor in much the same way as does the derivation of Zamolodchikov's c-theorem 
[6] which can be found for example in [7|. Zamolodchikov's c-theorem states that in a 1+1- 
dimensional unitary QFT the function c{r) defined by 



dc _ 3r3 
dr ~ 2~ 



((e(r)e(o))(M2)„ - (e)?K.) J , (1.5) 



with c(oo) = 0, has three fundamental properties: it is monotonically decreasing along RG-flows; 
it is positive or zero for all values of r and its derivative given by (jl.Sp vanishes only at fixed 
points r = r^:. At those points c{r:f) = c, the central charge of the CFT describing the critical 
point. These properties can be easily established from the definition (jl.Sp and the properties of 
the trace of the stress-energy tensor. Note that since r is an energy scale, the only fixed points 
are at r = and r = oo; but in certain cases, where the RG fiow approaches other fixed point 
at finite energy scales, a plateau structure is observed (see the famous example [8] for a related 
scaling function, and other works, e.g. [lOlIII]) where the c-function has been studied for various 
models). 

It is clear that for r = 0, oo the functions A(r) ()1.3p and c{r) (jl.Sp are proportional to 
each other through (jl.2p . However, given the rather different nature of the correlation functions 
involved one would naturally expect that 

A(r)^-^(n--] for r^0,oo. (1.6) 

24 \ nj 

Indeed for the Ising model, it is possible to obtain exact formulae both for c{r) and A(r) which 
confirm (|1.6p . Interestingly, they also show that the two functions are qualitatively very similar 
so that one may ask whether or not A(r) in general satisfies the properties of a c-function stated 
above. The first suggestion that a new c-theorem could be formulated for the function A(r) 
subject to negativity of (jl.3p appeared in the work [1]. Given the connection between the twist 
field T and the entropy, it is also interesting to mention the work [U] where a new c-function 
directly given in terms of the entanglement entropy was proposed. 

The idea that c{r) and A(r) should be qualitatively similar beyond the Ising example is 
strongly supported by results obtained in [12]. The similarity between the functions c{r) and 
A(r) can be best appreciated by comparing Fig. 1 in [10] to Fig. 1 in [12]. In particular, in |12j . 
a plateau structure is observed for A(r), with plateaus at the same energy scales as those of 
c{r). These results have provided the initial motivation to carry out the present investigation. 

The aim of this letter is to provide arguments supporting the conjecture that A(r) satisfies 
the properties of the c-function. 



2 Properties of A(r) 

Based on expected properties of unitary QFT, we will argue that 

(0(O>Mf - (e)K2 < V < r < oo. (2.1) 

Further, consider a model described by a CFT perturbed by a primary spinless field (p of con- 
formal dimensions A = A < 1/2, 

A = AcFT + gd'^z4>{z,z), (2.2) 

where Ac ft is the action of the underlying CFT on M?. Then, we will show that 

(0(0)^;? - (e)M2 oc r2^(^-i) (2.3) 

as r — ;• 0. 

From ()1.3p . Equation ()2.ip ensures that A(r) < 0, hence implies monotonicity of A(r). The 
property of positivity A(r) > for < r < cxd then automatically follows by integration of (jl.Sp . 



A(r) = --y^ dss((e(s))^n-(e)ffi2), (2.4) 

which converges by factorization of correlation functions at large distances and automatically 
implements A(oo) = 0. Equation ()2.ip says that A(r) can only vanish at r = or r = oo, the 
fixed points of the RG flow. It does so at r = oo thanks again to factorization at large distances. 
Thanks to (|2.3p and (|1.3p . it also does so at r = for all n, using the assumption that A < 1/2. 

It is interesting to note that these properties do not necessarily hold for operators other than 
the branch point twist field. The results of |1UJ (Figs. 2 and 3) and [TT] (Fig. 3) show explicit 
examples of fields for which the monotonicity property of A(r) does not hold. 

Zamolodchikov's positivity principle, used to prove the c-theorem, does not work here, so we 
develop new arguments: 1) analysis of the IR and UV regions using form factor expansions and 
perturbed CFT, respectively; 2) general QFT-based intuitive arguments. 

2.1 Neighbourhood of the IR and UV fixed points 

We begin by considering the large distance region. In the two-particle approximation we write 

-((e(r))^,.-(e)«.)-E /^^2^'"'(^^'^2)(^S7T^) e—^'"'''''^-'''\ 

(2.5) 
where for simplicity we assume that we have, on the right-hand side, n copies of a QFT with 
a single particle spectrum. The functions above are defined as Fg (^1,6*2) := {0\O{0)\0i92)ab-, 
where |0) is the vacuum, \6i02)ab is a two-particle asymptotic state, and 61^2 are rapidities. 
For integrable models, the twist field form factor is known to be 



sinh 



j^^^Mfi^lM) sinh ( M2(fc-a)-l)-g ^| F^ab^^iiTT, n) ' 



and the form factor. 



^e\ah^n n \ o 2_^^in(M) 



i^2 ' (^1,^2) = '^^r,{' rZnUM ^ (2-7) 



with 6 = 6i — 02 and {T) the vacuum expectation value of the twist field. The normalization of 
the form factor of G, F2 {itt) = 27rm'^, is fixed as explained in [T3|. Form factors with a ^ b 
are zero. F^^^{6,n) are the minimal form factors, analytic solutions to the equations 



F^?„(0, n) = Sab{e)F^i^{-e, n) = F^,^{27rin - 9, n), [21 



for any a,b = 1, . . . ,n. 

Inserting these expressions into ()2.5p . we can write 



m^sin- f Ko(2'mrcosh^) F^^ (xA) F^^ 



[x,n] 



27r 7 cosh --cos- F^.^{m,l) F^^^{i-K,nY 

—00 

where Ko^t) is a modified Bessel function. 

Clearly, the sign of ()2.9p is only determined by the minimal form factor product, as all the 
other quantities in the integrand are positive. It turns out that for integrable models of the type 
considered here, the minimal form factor always admits some integral representation of the form 



11 .. / dtf{t) ^., 



fOO 

Fj^|„(x, n) = exp / .^" V^' .. sin 
/o tsmh(nt) 



it I ix 

2 V vr 



(2.10) 



where fit) is a real function which depends on the scattering matrix of the model. The minimal 
form factor is in general a complex function, however the product 



11, X 1 1 . s ^ ; (it , , / 1 — cos — cosh nt 1 — cos — cosh t 



poo 

Fi!.(..i)Fi(.,„)-=expy^ i/w(^ - 'T^;;,; + ' ^2^; j. ("d 

is real and positive. This shows that near the infrared fixed point (mr large) the function A(r) 
defined in ()1.3p is negative. In addition, the presence of the exponential ()2.5p ensures that the 
value of the integral is larger for smaller values of mr. Note that for fields other than the branch 
point twist field, there is no reason to expect that the present argument, which depends on the 
particular form of the form factors, gives negativity. 

A model in which the correlation function involved in (jl.3p is known exactly is the Ising field 
theory. In this case, the two-particle approximation ()2.5p is exact, and we have [T5|[T6]: 

F^^^\e) = -2-Kim^smh-. (2.12) 

We may directly adapt our formula (|2.9p to the Ising case 



2 



00 



vr f i^o(2"ir cosh|)sinh^sinh| 



/,^/ Nv ,^v \ in" vr / , i\ol^"i^cosnTTismni5-smnTT 
n((e(r))^.-(e)M.)=-— ^o«^ / dx _, ^^1 _ f ^> (2-13) 



vr 2n / cosh ^ — cos r 



-00 



by employing F^^^{x,n) = — isinh^. This shows negativity for all < mr < cc. 

Let us now turn to the short distance behavior of (|1.4p . using (|2.2p . On dimensional grounds, 
the coupling constant g is related to a mass scale m as g ^ m?'~ , and we will take 5 > and 4> 
"positive" (see the next subsection) so that the spectrum of the theory is bounded from below. 
There is a direct relationship between the perturbing field and in the massive model, 

e{z, z) = 4vr5(l - /\)(t){z, z). (2.14) 



For A < 1 this equality is exact in the sense that no higher order corrections in g occur |17| . 

The expectation value (0('^))Aijf can be evaluated through the operator product expansion 
(OPE) of the left-hand side of (jl.4p . Keeping (j2.14p in mind, we write 

(0(r)r(O))(R2),^ _^^ {Of,)iR2)^ 

in terms of some fields O^ of the massive QFT. Considering the zeroth order of conformal 
perturbation theory [13], we directly replace the structure functions by their CFT value. The 
leading term of the expansion ()2.15p will involve a field Oq, written as the composite field : (pT :, 



((/.(r)r(0))(R2)„ ^ ^:^r:^2(A.,r:-A-Ar)iiiZlVk + . . . (2.16) 



where we have used dimensionality arguments to re-write C'?^'{r) = C^^ 'r^*- ■■*'''■ 'T' in 

terms of the conformal dimension of the new field A. ^7-., and a dimensionless constant C'fj- ' . 

It is possible to fix A^j^t": by comparing the OPE above to the standard CFT computation 
of a correlation function of the form: 

(4>{z,z)T{Q)0{x,x))(^2^ r2^(^-i) 



(2.17) 
where r = |z| = ^/z^ and f{z) = zn is the conformal transformation that unravels the Riemann 
sheets conformally mapping them to M^. The last equality is a short distance expansion. Here 
O is an arbitrary product of local fields, not necessarily primary, at positions represented by the 
sets x,x. From the comparison between ()2.16p and ()2.17p we can fix 

A:r0: = ^ + Ar and Cff- = ^. (2.18) 

We can then use dSlD, (PTTUj) and ^nE\i to rewrite the difference pT^i]) as 

2 ( aMmrf^^^-^'^ \ , , 

(0(O)ais - (©>R2 = "i' ^ J2A ^ + • • • (2.19) 



where ^'KgiX — A) = am?' ^ , (0)ir2 = [im?' and 



{T) 



/3m—, (2.20) 



and a, /3 and ^ are all dimensionless constants. With n > 1, this shows ()2.3p . 

Clearly a > if A < 1/2 and g > 0. Hence, negativity of (|2.19p at short distances requires 
/? < 0. Although expectation value {T)(w^\ is positive as it represents the partition function of 
the theory on the manifold A^q (the exact value for the Ising model is always positive [3]), we 
do not have a derivation of the negativity of (: (pT :)m2\ (it follows from the arguments of the 
next subsection). For the Ising model, (: (j)T :)(k2)^ may be evaluated explicitly using (j2.13p : 

n{{Q{r))M^i-{Q)R^)=-^cos^T{\-^\ {mr)^~' for mr « 1. (2.21) 

The leading mr behaviour is as expected (with A = i) and the overall sign is indeed negative. 



2.2 General arguments 



We have established that (j2.ip holds for large distances in a large class of integrable models, 
and we have shown (|2.3p in general. Further, formula (|2.13|) for the Ising model shows that (|2.1|) 
holds for any mr in this case. We now provide model-independent arguments, based on expected 
physical properties of unitary models, strongly suggesting that (|2.1|) holds for arbitrary values 
of mr. Here, we use in an essential way the geometric interpretation of the branch point twist 
field, hence these arguments do not apply to any other field. 

Note that proving (|2.ip (for n > 1) is equivalent to showing that 

^ mr))M-) < 0. (2.22) 

Indeed, if (©(?")) Aljf decreases with n, then so does (©(r))^^" — (0)r2. Hence, the largest value 
of the latter for n > 1 is at n = 1, where it is because A^q^^ = M^. Similarly, if 

^ ((0(O)^S) > 0, (2.23) 

then (|2.ip follows, because of factorization of correlation functions at large r. In particular, 
establishing ()2.23p immediately shows negativity of the coefficient /3 in ()2.20p . 

Our main argument uses the idea of virtual particle propagation. We re-interpret unitarity 
as "positivity" of the perturbing field (j) (hence of 0): (p should be an appropriate normal-ordered 
(i.e. renormalized) product of an operator tp and its hermitian conjugate, (/> = (ip^'ip), in analogy 
with the factorization of positive-definite matrices. Then contributions to the expectation value 
(0('"))ai" come from virtual particles created and annihilated at the point (r, 0), and propagating 
in A4q- Every path contributes a positive amplitude proportional to the exponential of minus the 
single-particle Euclidean action (i.e. the Brownian motion measure of the path), with possible 
branching due to interactions. Some of these paths go around the origin. As the angle around 
the origin 27rn increases, these self-interaction contributions become less important, because 
the distance traveled is greater. Whence the derivative with respect to n is negative, giving 
()2.22p . A similar argument leads to ()2.23p for n > 1. Indeed, as r decreases, more and more 
self-interaction loops must travel around the origin, hence giving lesser contributions. 

A way to study the self-interaction loops around the origin is to use angular quantization. 
Let us consider as an example the Klein-Gordon theory, and explicitly show ()2.22p in this case. 
Angular quantization was developed quite generally in |201 \2T\ in the context of form factors in 
integrable models; the Klein-Gordon angular quantization described in [21] allows us to evaluate 
correlation functions. The construction of the branch point twist fields in angular quantization 
was described in [3]. Let us summarize few key ingredients. 

We are interested in the operator G oc : (^^ : where ip is the Klein-Gordon field; the normal- 
ordering is a point-splitting regularization, with a subtraction proportional to the identity. We 
first compute the two-point function (c^(r, 0)c^(r',0)), then take the limit r — )• r' . In the angu- 
lar quantization approach, correlation functions are expressed as traces over the space of field 
configurations on the half-line (representation denoted by vr^). The density matrix used in this 
trace is the operator performing a rotation by the angle necessary to go around the origin. The 
corresponding conserved charge associated to rotation, denoted by K, is the Hamiltonian of the 
theory. The presence of the branch-point twist field means that the angle around the origin is 
2?™. Hence, the density matrix is e^'^*"^: 

, , _ Tr.,(e^-^vr^(---)) 

{■■■)M^S- Tr^ (e27rin.A-) • (2.24) 



For the Klein-Gordon theory, ttz is a representation of the Heisenberg algebra with oscilla- 
tors hi, that satisfy [6;/,6i/'] = 2sinh(7rz^)(5(i^ + v'). Then, the following relation holds [3]: 

"^ T>r (e^'^'"-^^) ^ ^'^'^'^ s\nh{^nl) ^^'^ + ^')- Further, the bosonic field is expressed as is [21] 

7rz(^(r,0)) = ^ / duKKUmr) with K^mr) = - d^e"™^ ^"^'^ V'^^ (2.25) 

Employing all these definitions, the two-point function can be written as 

(^(r,0)<^(r',0))^n = - r due-''''^^^^^K,,{mr)K.Umr'). (2.26) 

^ J-oo sinh(7rni/j 

Since the conformal point is a free boson, this function diverges logarithmically when r ^ r' . 
However, this divergence is independent of n, whence we differentiate then take r = r': 

'^'- ""^'f '-^^^ = -4 r d.. "t^""\ |K..(.^r)|^ < for ,nr / 0. (2.27) 

an J -OO sinh (vrni/) 

This establishes (|2.22|) for the Klein-Gordon theory. 

For more general, unitary models, the argument goes as follows. The operator K, as it does 
above, should have positive imaginary eigenvalues in order for the trace (j2.24p to be well-defined. 
Hence, let us write iK = —J for a positive operator J. Differentiating with respect to n, we find 

— (</)(r,0))^n = -27r((J0(r,O))A4;f - (^)^j(</'(r-,0))A^jj) . (2.28) 

Since the measure is rotation invariant, and since J is proportional to the generator of rotations, 
we have ([J, (/)(r, 0)]);v(jf = 0. Hence, J and i;^(r-, 0) can be interpreted as "classical" statistical 
variables, and the derivative with respect to n is the negative of their statistical correlation. 
We expect this statistical correlation to be positive: the statistical variable J is an "energy", 
composed of a kinetic energy (the conformal part) and a potential energy V (the perturbation 
by (j)). Indeed, a moment's thought shows that if the average potential energy {V)j at fixed 
total energy J increases with J, as should be expected, then J is positively correlated with V . 
Finally, let us further justify the latter angular-quantization argument through a drastic 
simplification. Instead of propagating a half-line around the origin for an angle of 27rn, we 
reduce to a finite number of degrees of freedom: we consider the propagation of a quantum 
mechanical particle along a circle of circumference 2im. This simplification is expected to 
provide the right sign of the variation with respect to n, which comes from particles propagating 
around the origin. The operator J is replaced by the Hamiltonian H of the quantum system, 
and the perturbing field cj){r) is replaced by the potential energy V . The trace becomes 

j dx {x\e-'^-"'^^\x) ' ^ ' ^ 

Since quantum mechanics in imaginary time corresponds to a stochastic problem, we need to 
evaluate the average of the potential V{x), with an un-normalized measure given, for any value 
of the position x, by the probability for a random walk in that potential to start and end at 
rr in a time 27m (more precisely, to come back to a position in a small, fixed neighborhood 
[x — 6,x + 6]). As time increases, this probability decreases for any x. However, at lower values 
of the potential, nearer to the absolute minimum, the additional time given to the particle is 
more likely be spent near to its original position than it is at larger values, because the particle 
has a tendency to fall back to the minimum of the potential. Hence, as time increases, lower 
values of the potential get more relative weights. This implies that the average of the potential 
decreases as the time 27rn increases. 



3 Conclusion 

24nA(r) 

In this paper we have provided evidence that the function c(r) = — ^ with A(r) given by 

^_^ n^ — 1 

(j2.4|) satisfies Zamolodchikov's c-theorem even though it is not Zamolodchikov's c-function. 

We employed several types of arguments and techniques to establish our main conclusion. 
First, we used a form factor expansion to show that A(r) is monotonically decreasing for large 
distances; second we used the OPE of operators and T to investigate the short distance behav- 
ior of A(r). For the Ising and Klein-Gordon models we proved that A(r) is also monotonically 
decreasing for all r. We then argued, from various physical ideas, that this holds for general 
unitary models. It would be very interesting to put these arguments on a more solid basis by, 
for instance, explicit perturbative calculations. 
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